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In this work, we have studied (anti)neutrino induced charged current quasielastic scattering from 
some nuclear targets in the energy region of E,, < 1 GeV. Our aim is to confront electron and 
muon production cross sections relevant for -o- Uf, or o v,, oscillation experiments. The effects 
due to lepton mass and its kinematic implications, radiative corrections, second class currents and 
uncertainties in the axial and pseudoscalar form factors are calculated for (anti)neutrino induced 
reaction cross sections on free nucleon as well as the nucleons bound in a nucleus where nuclear 
medium effects influence the cross section. For the nuclear medium effects we have taken some 
versions of Fermi gas model(FGM) available in literature. The results for (anti)neutrino-nucleus 
scattering cross section per interacting nucleons are compared with the corresponding results in free 
nucleon case. 

PACS numbers: 12.15.Lk,12.15.-y,13.15+g,13.60Rj,21.60.Jz,24.10Cn,25.30Pt 


I. INTRODUCTION 

With the measurement of 0i3 at the nuclear reactors [H-Q and finding a clear evidence for Pe disappearance, now 
the emphasis is upon determining mass hierarchy in the neutrino sector as well as to find simals of CP violation in the 
leptonic sector Si- For these physics goals, accelerator experiments like T2K 0, NOz/A etc. are taking data and 
looking for o z/g or o zzg oscillation signals and experiments like LBNO Q, DUNE T2HK [iflllll) etc. are 
planned. Most of these oscillation experiments are being performed in the (anti)neutrino energy region of ^ 1 GeV. 
In this energy region the major contribution to the event rates comes from the charged current quasielastic(CCQE) 
lepton production process followed by charged current induced one pion production process. There are other channels 
which may contribute to the lepton event rates like strange particle production through |A5'|=1 processes but they 
are Cabibbo suppressed, while their production through strangeness conserving |A5'|=0 processes are suppressed due 
to the threshold effect. The contribution to the lepton event rates from the deep inelastic scattering is expected to 
be small. The importance of CCQE lepton production process is two fold. Eirstly, the observation of charged lepton 
in the final state is the cleanest signature of a (anti)neutrino interaction and secondly this is the simplest process 
using which energies of the incoming (anti)neutrinos may be determined. However, there are theoretical uncertainties 
involved while inferring the energy of neutrinos even from CCQE processes on nuclear targets and this has been 
discussed in detail recently by many authors [l^ - [T^ . 

In the electron and muon (anti)neutrinos processes, the contribution to the cross section from CCQE process would 
be different for the production of electrons and muons, due to lepton mass and other effects even in the presence of 
the universality of weak interaction. Recently, Day and McEarland studied the effect of lepton mass, radiative 
corrections and uncertainties in the nucleon electroweak form factors including the second class currents(SCC) on the 
(anti)neutrino CCQE scattering cross sections from the nucleon targets. In their work, it has been shown that the 
radiative corrections at the tree level CCQE process may lead to important difference between electron and muon 
production cross sections [l^, as it is proportional to log('^), where E* is the outgoing lepton energy in the center 
of mass frame and mi is the mass of the charged lepton. Furthermore, the variation in the axial dipole mass Ma 
which has been recently discussed in literature will also lead to difference in electron and muon cross sections. The 
variation in the experimental measurements of Ma has been recently found out to be quite large i.e. 0.99 GeV to 
1.35 GeV [2ll - l^ . from the world average(1.026 ± 0.021 GeV 0 and 1.014 ± 0.014 GeV [s^). A change of 10% in 
the value of Ma say from 1 GeV to 1.1 GeV would result in almost an increase of 10% in the CCQE scattering cross 
section. Generally, the pseudoscalar form factor Fp{Q^) is expressed in terms of axial vector form factor Fa{Q^) 
using the Goldberger-Trieman relation and PCAC. However, other parameterizations of the pseudoscalar form factor 
Fp{Q‘^) based on Chiral Perturbation Theory and Lattice calculations have also been discussed in literature 
Moreover, if one considers the presence of second class currents, then there would be additional contribution to the 
(anti)neutrino nucleon cross sections due to the form factors F^{Q^), Fj^{Q‘^) associated with them. The inclusion or 
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noninclusion of these contributions may translate into the systematic uncertainties in the determination of event rates. 
Some of these effects have not been taken into account in the (anti)neutrino Monte Carlo generators like GENIE [s^, 
NEUT NUANCE [i^, NuWro GiBUU HI, etc. 

In almost all the present generation (anti)neutrino experiments, moderate to heavy nuclear targets like 

etc. are being used. Recently, the measurements have been performed for distribution [2ll - [^ . 
double-differential scattering cross sections and total scattering cross sections [12, [13, HE 113, using some 

of these nuclear targets. In these nuclear targets nuclear medium effects play an important role and the effect has been 
found to be substantial in the low energy region. Most of the present generation Monte Carlo generators are using 
relativistic Fermi gas model given by Smith and Moniz (43j| or the model discussed by Llewellyn Smith H . However, 
there are other variants of Fermi gas model available in literature like that of Gaisser and O’Connell [^, Singh and 
Oset [i^,|43|, Nieves et al. [iEl, etc. which have been applied to study (anti)neutrino-nucleus reactions. In addition to 
these works many other theoretical models like superscaling approach [49j . mean field approximation [HE - relativistic 
meson-nucleon model [5l|, relativistic Green’s function model plane wave impulse approximation(PWIA) [HEHE, 

distorted wave impulse approximation(DWIA) [HE, etc., have been used in literature. Moreover, in the context of 
Fermi gas models the effect of two particle-two hole(2p-2h) correlation, meson exchange currents and multinucleon 
mechanism have also been considered in literature [56l - l59l| . These studies have been recently summarized in some 
review articles [lE [HE [Ml , but they have not been incorporated into the present Monte Carlo generators being used 
in the neutrino oscillation experiments [63 |. 

In the present work, we have studied vi, vi] {I = e, induced CCQE scattering from some nuclear targets like 

and in the energy region of < 1 GeV including the effect due to lepton mass and its kinematic 

implications, radiative corrections, form factors, second class currents, etc. We have performed the calculations using 
Local Fermi Gas Model(LFG) [HE - lbE and also compared the numerical results with the different Fermi gas models of 
Smith and Moniz [iE, Llewellyn Smith [iE, Gaisser and O’Connell [IE. Using the different nuclear models we have 
studied the difference in the lepton(electron vs muon) cross sections due to the axial dipole mass, pseudoscalar form 
factor, radiative corrections and effect of second class currents. Furthermore, in the Local Fermi Gas Model [gE - IHE 
we have also included the nucleon-nucleon interactions due to which response of electroweak transition strength is 
modified. These modifications are calculated by incorporatin g th e interactions of particle-hole(lp-lh) excitation in 
the nuclear medium in a random phase approximation(RPA) [dE IHE ■ Our aim of this work is to study the role of 
dominant nuclear medium effects like Fermi motion, binding energy and nucleon correlations on the various physics 
inputs which lead to difference in electron and muon scattering cross sections in the case of per nucleon target as 
discussed by Day and McFarland [IE. In view of this, the effect of 2p-2h, meson exchange currents and multinucleon 
mechanism have not been considered. This is at present beyond the scope of this paper and it may be studied as a 
separate work. In section-HIl we present the formalism in brief, in section- HlIl results and discussions are presented 
and the findings are summarized in section- HVl 


II. FORMALISM 

The basic reaction for the quasielastic process is a (anti)neutrino interacting with a (proton)neutron target given 

by 


vi{k) + n(p) ^ I {k') + p{p') 1 , 

vi{k) + p{p) l^{k') + n{p') J 


( 1 ) 


where k, k' are the four momenta of incoming and outgoing lepton and p, p' are the four momenta of initial and 
final nucleon, respectively. The invariant matrix element for the charged current reaction of (anti)neutrino, given by 
Eq.([T]) is written as 


M = cos 9c L 


( 2 ) 


where Gp is the Fermi coupling constant (=1.16639x10 ^ GeV ^), 9c{= 13.1°) is the Gabibbo angle. The leptonic 
weak current is given by 


If, = u{k')^f,{l±^^)u{k), 

where (-|-ve)—ve sign is for (antineutrino)neutrino. is the hadronic current given by 


( 3 ) 


= u{p')T^u{p), 


( 4 ) 
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with 


= FY{Q^)r + F^{Q^)za^''j^+F^{Q^)^ 

+ Fa(Q^)i^1^ + Fpm^j^ + F3^(Q2) (P±|T^5^ 


( 5 ) 


Q^{= —q^) > 0 is the four momentum transfer square and M is the nucleon mass. F^ 2 iQ^) the isovector vector 
form factors and Fa{Q^), Fp{Q^) are the axial and pseudoscalar form factors, respectively. F^(Q^) and Fj^{Q‘^) 
are the form factors related with second class current. Using the leptonic and hadronic currents given in Eq.(j31) and 
Eq.Q, the matrix element square is obtained by using Eq.@ as 


1-^1^ = ^cos^^c 


is the leptonic tensor calculated to be 




where 


/J.lv' 


= 8 k ■ fc'] and 


= 8 k'^k^, 

where the + sign(— sign) is for neutrino (antineutrino). 

The hadronic tensor is given by: 


( 6 ) 

(7) 

( 8 ) 


= ^Tr 
2 


{i)' + M)T^^{i) + M)f'- 


where = 7 O 

The hadronic current contains isovector vector form factors FY 2 {Q‘^) of the nucleons, which are given as 

FYaiQ^) = Fl^m - Fl^iQ^) 


(9) 


( 10 ) 


where F^^^\q'^) and F 2 ^^\q‘^) are the Dirac and Pauli form factors of proton(neutron) which in turn are expressed 
in terms of the experimentally determined Sach’s electric G^"(Q^) and magnetic G^{Q'^) form factors as 


F^iQ^) = 1 + 


F|’"(Q 2 ) = 1 + 


4M2 

4M2 


-1 


Gr(Q") + 




4M2 

[G"Af(Q^)-Gr(Q")] 


G^”(Q^) 


( 11 ) 

(12) 


G^"(( 5 ^) and G^{Q'^) are taken from different parameterizations [67 

The isovector axial form factor is obtained from the quasielastic neutrino and antineutrino scattering as well as 
from pion electroproduction data and is parameterized as 


Fa{Q^)=Fa{0) 


1 




Ml, 


-2 


; F^(0) =-1.267. 


(13) 


The pseudoscalar form factor is determined by using PCAC which gives a relation between Fp{Q^) and pion-nucleon 
form factor gjrNNiQ'^) and is given by [43| : 


ml g-nNNiQ'^) 


2 . ^ 2M^Fa{ 0) [FAm _ 

{ Fa{ 0) (ml + Q^) g^NNiO) 


(14) 


where is the pion mass and g-^NNiO) is the pion-nucleon strong coupling constant. 

Fp{Q^) is dominated by the pion pole and is given in terms of axial vector form factor Fa{Q^) using the Goldberger- 
Treiman(GT) relation 


Fpm = 


2M^Fa{Q^) 

to 2 _|_ Q 2 


(15) 
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The form of pseudoscalar form factor Fp{Q‘^) using PCAC may also be written as 


fM) 


M W 2mlF^ \ 
Q2 [\rnl + Q^) 


(MFa{Q) 

V 


ff7i-Ar7v(0)A(5^ \ 

to 2 j 


+ 2 MFa(Q 2 ) , 


where = 13.21, F^ = 92.42 MeV and A = 1 + 

Pseudoscalar form factor using Chiral Perturbation Theory(CliPT) is given by 


ml + Q^ +-3- 


(16) 


(17) 


where axial radius ta = 

The form factors F^(Q^) and Fj^{Q'^) are associated with the second class current(SCC). There are no compelling 
reasons for their existence as they violate charge or time symmetry and in the case of F^{Q^) also the conserved 
vector current. Almost all the current calculations of neutrino reactions assume SCC to be zero. However, there are 
some experimental analyses of semileptonic weak interactions like beta decays, muon capture and neutrino scattering 
in the lA^j = 0 sector which give u ppe r limits on these form factors which are consistent with the constraints of the 
present data on these processes l72l - l^ . In Ref. , the upper limit for second class vector current obtained from 
neutrino experiments is 1.9. We have used the following expressions for F^ {Q^) as given in Ref. 721 


F^{Q^) = 


F^ 


( 0 ) 




(18) 


To observe the maximum effect of second class vector current we have taken ^3^(0) = 1.6 on the upper side of the 
limit with = 1 GeV for our numerical calculations. Another expression for F^ (Q^) as given in Ref.[2^ is 

F^{Q^)=4AF^iQ^). (19) 

The axial form factor associated with the second class current F^{Q‘^) is taken as [20l [t^ 

F3^(Q2) = 0.15 FAiQ^). (20) 


The parameterization of form factors discussed above will be used in the evaluation of the CCQE cross section. The 
differential scattering cross section for reaction given in Eq. o in the laboratory frame is in general written as, 


da = 


{2TT)‘^6^{k + p — p'— k') d?k' 


d?p' 


-EE|7Wp. 


4(fc-p) {2tt)^2Ei {2tt)^2Ep 

The double differential cross section afree{Ei, fl;) on free nucleon is then obtained as 


( 21 ) 


^ free{Fp 11 /) 


d^a 

dEi didi 


w\ 

Qdn'^E^EnEp 


SE|7W|2(5[go + A„ 


Ep] 


( 22 ) 


Inside the nucleus, the neutrino scatters from a neutron moving in a finite nucleus of neutron density p„(r), with 
a local occupation number n„(p, of the initial nucleon of momentum p localized in a region r, the radius of the 
nucleus. In the local density approximation the scattering cross section is written as 


a{Ei,ni) 


2drdp-—-^nn[p, r)af 

ree (E/,1!/) 


(23) 


where afreeiEi,fli) is given by Ea. (l22l) . The neutron energy En and proton energy Ep are replaced by E„(|p|) and 
Epi\p + qI), where p is now the momentum of the target nucleon inside the nucleus. This is because inside the 
nucleus the nucleons are not free and their momenta are constrained to satisfy the Pauli principle, i.e., p < pp^ and 
p'{= \p + g\) > PFp, where pp^ and pp^ are the local Fermi momenta of neutrons and protons at the interaction point 

in the nucleus and are given by pp^ = [37r2p„(r)] ^ and pp^ = [37r^/9p(r)] p„(r) and Pp{r) are the neutron and 

proton nuclear densities which are given in terms of the nuclear density p{r) : 


Pn{r) 


Pp{r) = ^P{r). 


(24) 
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Cl C 2 Q—value B.E. pp 


Cl C^ C2 u u V V 



1.692 

1.692 

1.082* 

1.082* 

16.8 

13.9 

25 

221 

221 

^°Ar 

3.64 

3.47 

0.569 

0.569 

2.5 

8.0 

30 

259 

242 


4.05 

3.971 

0.5935 

0.5935 

6.8 

4.8 

36 

263 

251 

20Spb 

6.89 

6.624 

0.549 

0.549 

2.4 

5.5 

44 

283 

245 


TABLE I: Different parameters used for numerical calculations for various nuclei. Ci and C 2 are the density parameters(in 
Fermi units) defined for modified harmonic oscillator as p(r) = po(l + C 2 (^)^)ea:p( —(^)^) and for 2-parameter Fermi density 
as p(r) = po/(l+ea;p( ^~^^^ ))■ For we have used modified harmonic oscillator density(* ci is dimensionless) and for *°Ar,^^Fe 
and nuclei, 2-parameter Fermi density have been used, where superscript n and p in density parameters(c"’^; i=l,2) 

stand for neutron and proton, respectively. The Q—value of the reaction, binding energy {B.E) and Fermi momentum(pF) for 
different nuclei are given in MeV. 




FIG. 1: Diagrammatic representation of the neutrino self-energy corresponding to the ph-excitation leading to vi+n^ l~ + p 
in nuclei. In the large mass limit of the intermediate vector boson(i.e. Mw —t oo) the diagram (a) is reduced to (b) which is 
used to calculate |Af P in Eq.®. 


The density parameters have been taken from Ref. [zllii and are summarized in Table-lH For the antineutrino 
induced reaction on free nucleon or nucleons bound in a nucleus the role of neutron and proton get interchanged. 
Furthermore, in nuclei the threshold value of the reaction i.e. the Q—value of the reaction(Qr) has to be taken into 
account, which we have taken to be the value corresponding to the lowest allowed Fermi transition or Gammow-Teller 
transition. 

These considerations lead to a modification in the 6 function used in Ea. (l2^ i.e. i5[(?o + En — Ep] is modified to 
i5[(7o + En{p) — Ep{p + q) — Q^] and the factor 

occurring in Ea. (E^ is replaced by — (l/7r)Iml7jv(go, where UN{qQ,q) is the Lindhard function corresponding to 

the particle hole(ph) excitation shown in Fig.dT]) and is given by: 


UN{qo,^ 


dp rinip) [1 - np{p + ^] 

(27r)3 EnEp qo + En{p) - Ep{p + q) + ie 


( 26 ) 


where qQ=E^ — Ei — Qr- For the antineutrino reaction the sufhx n and p will get interchanged. 
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The imaginary part of the Lindhard function is obtained to be [47 


ImllNiqo,^ 


1 M 2 
2tt 


[Ef, - A] 


(27) 


with > 0, Ef 2 — qo < Epi and 


-go+|giyi+^ 


< Epi, where Ep^ = sJvFn^ + ^ Ef 2 = \/pFj,'^ + M^ and 


A = Max 


M, Ef 2 — qo, 


-<zo + |g1^/l + ^ 


With inclusion of these nuclear effects, the total cross section a{E^) is written as 


a{E^) = 


X 


rrmc 

.X l-k'm. 


2Gf^ cos^ 0c I 

r^dr 1 

k'dk' / 


■Jk'mir 

-'QLn 


L^,J>^MmUN[E, - El-Qr,^. 



(28) 


In the above expression Tmin and Xmax are the minimum and maximum limits of nuclear size. In principle r should 
vary from 0 to oo but in the local density approximation density dies out at around 10 /m, therefore, for our numerical 
calculations we have taken integration limits for nuclear size from 0 to 10 fm. and k'^^^ are minimum and 

maximum values of outgoing lepton momenta. The energy and momentum of the outgoing lepton get modified due to 
the Coulomb interaction, which is taken into account in an modified effective momentum approximation(MEMA) [T7I |. 

In the local density approximation, the effective energy of the lepton in the Coulomb field of the final nucleus is 
given by ISItII: 


Eeff = El + Vc(r), 


where 


Vc{r) = Ma Zf J 




Z 


/ 


Ppir'] 


Z 


r' dr' 


f 


(29) 


with a as fine structure constant and Zf as the charge of outgoing lepton, taken as —1 for neutrino and +1 for 
antineutrino. This leads to a change in the imaginary part of the Lindhard function occurring in Eq. 


ImUN[Ei, - El - Qr,^ -t ImUNiEv - Ei - Qr - Vcir),cf) 


Furthermore, in a nucleus the response of electroweak strength may change due to the presence of strongly interacting 
nucleons. These changes are calculated by considering the interaction of ph excitations in the nuclear medium in 
Random Phase Approximation (RPA) as shown in Figl^l The diagram shown in Figl2] simulates the effects of the 
strongly interacting nuclear medium at the weak vertex. The ph-ph interaction is shown by the wavy line in Fig|^ 
and is described by the n and p exchanges modulated by the effect of short range correlations. The effect of the A 
degrees of freedom in the nuclear medium is included in the calculation of the RPA response by considering the effect 
of ph-Ah and Ah-Ah excitations as shown in Figl^Kb). This is done by replacing Un by Un = Un + Ua, where Ua 
is the Lindhard function for Ah excitation in the nuclear medium and the expressions for Un and Ua are taken from 
Ref. 1^. The different couplings of N and A are incorporated in Un and Ua- The details of which are discussed 
in Refs, id^ . In the Appendix we have given the expression for the hadronic tensor in covariant form as well as 
the expression for the hadronic tensor when RPA corrections are incorporated. We must point out that the RPA 
corrections are implemented in the leading order terms only. 

Thus, in a local density approximation in the presence of nuclear medium effects including the RPA effect, the total 
cross section (j{Ey), is written as 


a{Ei,) = -2Gf^ cos^ 6c [ r'^dr f k'dk' f dQ"^ ^ 

Jr-min Jk'min. Qmin 

X ImUN[E^-Ei-Qr-Vc{r),^ (30) 

where is the modified hadronic tensor when RPA effect is incorporated and the energy transferred to the 

hadronic tensor also gets modified from q^ = E^ — Ei to qo = E^ — Ei — — Vc- Explicit expression of is given 

in the Appendix. 
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FIG. 2: Many body Feynman diagrams (drawn in the limit Mw —t oo) accounting for the medium polarization effects 
contributing to the process r'; + n —>■ + p. 


Now we discuss in brief the form of the various other Fermi gas model used in literature IdHj ] . In the Llewellyn 
Smith Fermi gas model , the cross section per nucleon in a nucleus is equal to the cross section for a free nucleon 
i.e. a free defined in Eql22l multiplied by [l — > where 


D 


Z for 2x < u — V 


0 for X > u + V 




for u — v<x<u + v 


(31) 


with X = y = (2^)i/3^ y = (22^1/3 N{= A — Z), Z, A are neutron, proton and mass numbers of the initial 

nucleus, respectively . t>f is the Fermi momentum and the three momentum transfer = \/■ 

Smith and Moniz [43| used the following expression for the double differential cross section in the Fermi gas model: 


(Pa 

dk'dVLi 


G'^k''^ cos^(i 


27r2M 


( 2 X) [ 


ITl? 1 . 2 A 


VF2 + [2lFi + ^IF„]tan^(-x) 


+ {Wp + Ws)mf /{MEl cos^(-x)) - 2WslM tan(-x) 

X sec(ix)[-Q^ cos^(^x) + l^^ sin 2 (ix)+wf]^ 


(32) 


where cosx = cos 9. The form of Wfs and other details are given in Ref. (d^ . 

Gaisser and O’Connell [d^ have used relativistic response function R{q, go), in a Fermi gas model to take into 
account nuclear medium effects, the expression for the double differential scattering cross section is given by 


dfa 

dilidEi 

R{q, qo) 


= C 


da fr 


d^i 


1 


Riq,qo), 

" d^pN 




EnEn' 
X e{\pN + pf^,), 


S{En + qo 


Eb 


En')9{pfn 


ItvI) 


(33) 


where is the Fermi momentum for the initial nucleon, N, iV'=n or p and C = A— Z for neutrino induced process 
and C = Z ior the antineutrino induced process. differential scattering cross section for the (anti)neutrino 

reaction on free (proton)neutron target and we have used the same expression for the form factors as used in the 
LFG for the numerical calculations. Different parameters associated with nuclear densities, Q^value of the reaction, 
binding energy and Fermi momentum used in the numerical calculations are summarized in Table-lH 















FIG. 3: Ratio vs E^, for neutrino(Left panel) and antineutrino(Right panel) induced processes in and 

The solid(dashed) line represent cross section obtained from electron(muon) type neutrino and antineutrino beams. For 
neutrino induced process N = A — Z, m neutron number and for antineutrino induced process = Z, is proton number, a a is 
cross section in nuclear target and has been evaluated using Local Fermi Gas Model(LFG) and LFG with RPA effect(LFG+RPA) 
and (j/ree is the cross section for the free nucleon case. 


III. RESULTS AND DISCUSSION 
A. Nuclear model dependence 

In this section, we present the results and discuss the findings. In Fig. |31 the results are presented for the ratio of 
scattering cross section per interacting nucleon to the scattering cross section on free nucleon target for (anti)neutrino 
induced processes in and in the energy region from threshold to 0.8 GeV. The results are 

obtained using Local Fermi Gas Model(LFG) i.e. the expression given in Eql28] and the Local Fermi Gas Model 
with RPA effect (LFG+RPA) i.e. using EplSHlto the cross section obtained for the free nucleon case using Eg l?I] on 
neutron(proton) target induced by neutrino (antineutrino). A similar study of v{v) — A cross sections in the present 
model at low energies(< 0.5 GeV) have also been performed by Kosmas and Oset in several nuclear targets 
including '^^Ar. The results obtained in this paper are in agreement with their work. 

Performing calculations using LFG, we find that in the nuclear medium effects like Fermi motion, Pauli 
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FIG. 4: The fractional suppression in cross section (5(JModei(= - Model ^ where a free is the cross section obtained for 

^ free 

free nucleon and OModei is per interacting nucleon cross section in Ar obtained by using different nuclear models. The results 
are presented for the cross sections obtained from different models of Fermi gas{aModei) viz. Smith and Moniz (dashed 

dotted line), Llewellyn Smith (dashed line), Gaisser O’ Connell [4^ (solid line), and with(double dashed dotted line) & 
without RPA(dashed double dotted line) effect using Local Fermi Gas Model. The top panel is for neutrino and bottom panel 
is for antineutrino induced processes. 


blocking, binding energy, result in the reduction of cross section by ~ 30(42)% at E^, = 0.3 GeV and around 20(30)% 
at = 0.6 GeV from free nncleon case for r'e(zze) induced processes. Inclusion of RPA correlation in LFG, reduces the 
cross section for Ve^Ve) scattering from free nucleon by ^ 55(56)% at E^ = 0.3 GeV and 35(45)% at E^ = 0.6 GeV. 
Similar results may be observed for ^^Fe and nuclear targets. In general, the reduction in the cross 

section increases with the increase in mass number. For and induced processes at lower energies the reduction 
is larger and for E^ > 0.4 GeV, the reduction in r'e(Pe) and cross sections is almost the same. This will be 

discussed separately in the next section when we compare electron and muon scattering cross sections. 

To compare our results with other variants of Fermi gas model, we have obtained total scattering cross section 
in ‘^^Ar using Fermi gas model of Smith and Moniz [i^, Llewellyn Smith and Gaisser and O’Gonnell [l^ and 
calculated fractional difference SaModei(= ), the results for which are shown in FigHl Here afree stands 

' ^free ' •' 

for the (anti)neutrino induced interaction cross section on free nucleon target and aModei stands for the (anti)neutrino 
induced interaction cross section for the nucleons bound inside the nucleus. The results for neutrino is different from 
antineutrino and is mainly due to the interference terms with Fa which come with an opposite sign. In the case of 
LFG with RPA effects, the effect of renormalization is large and this suppresses the terms with F 2 and Fa which 
results in a large change in neutrino vs antineutrino results. We find appreciable difference in the results when various 
nuclear models are used. 

For example, when calculations are performed by using Fermi gas model of Llewellyn Smith [dj, the cross sections 
get reduced from the free nucleon case by ~ 16(45)% at E^, = 0.3 GeV and around 8(30)% at = 0.6 GeV for 
Veipe) induced scattering processes. While when one uses Fermi gas model of Smith and Moniz [43l| this reduction 
in the cross section from the free nucleon case is ~ 42(45)% at E^, = 0.3 GeV and ~ 24(22)% at Ei, = 0.6 GeV. 
When the calculations are performed using Fermi gas model of Gaisser and O’ Connell [d^ the reduction from the free 
nucleon case is ^ 38(52)% at E,, = 0.3 GeV which becomes 20(35)% at E„ = 0.6 GeV for Vepe) induced processes. 
Performing calculations in the LFG, reduces the cross sections by ^ 25(44)% at E„ = 0.3 GeV and around 15(30)% 
at E,, = 0.6 GeV for Vepe) induced processes. Including RPA correlation with LFG, reduces the cross section by 
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FIG. 5: Ai = £qj. neutrinofleft panel) and antinentrino(right panel) induced processes in and Ar targets. 

’^ueiv'e.) 

Here / stands for the resnlts of the cross sections obtained (i) for the free nncleon case(solid line) (ii) in the Local Fermi Gas 
Model(dashed line) and (iii) for LFG with RPA effect(dashed dotted line). 


~ 54(58)% at El, = 0.3 GeV and 32(45)% at Ey = 0.6 GeV. The nuclear model dependence is found to be larger in 
the case of scattering than in the case of Ve{ve) scattering in the energy region of Ei, < 0.8 GeV. 

For induced scattering processes, reduction in the cross section from the free nucleon case is ^ 18(47)% 

at iJy = 0.3 GeV and ^ 7(30)% at Ey = 0.6 GeV, when cross sections are obtained using Fermi gas model of 
Llewellyn Smith . Fermi gas model of Smith and Moniz reduces the cross section from the free nucleon case by 
~ 40(30)% at Ey = 0.3 GeV and around 22(12)% at Ey = 0.6 GeV. When calculations are performed by using Fermi 
gas model of Gaisser and O’ Connell [i^ the reduction from the free nucleon case is ^ 40(52)% at Ey = 0.3 GeV 
which becomes 20(35)% at Ey = 0.6 GeV. The calculations are also performed in the LFG and we find that the 
cross sections get reduced by ^ 26(45)% at Ey = 0.3 GeV and around 15(30)% at Ey = 0.6 GeV for induced 

processes. Inclusion of RPA correlation with LFG, reduces the cross section for scattering from free nucleon value 
by ^ 58% at Ey = 0.3 GeV and 32% at Ey = 0.6 GeV, which for induced scattering is ^ 60% at Ey = 0.3 GeV 
and ^ 45% at Ey = 0.6 GeV. 


B. Effect of lepton mass and its kinematic implications 

There are two types of corrections which appear when lepton mass mi{l = e,/i) is taken into account in the cross 

section calculations for the reaction i^;(7';) + iV —>■ + N', {N,N' = n,p) which can be classified as kinematical 

and dynamical in origin. The kinematical effects arise due to Ei ^ |fc'| in presence of uii and the minimum and 

maximum values of four momentum transfer square {Q^ = —q^ > 0) i.e. and Q^ax g®ts modified, affecting 

the calculations of total cross sections. These effects are negligible for highly relativistic leptons but could become 

important at low energies near threshold specially for muons. On the other hand, the dynamical corrections arise as 

2 

additional terms proportional to in the existing contribution of vector and axial vector form factors as well as 
new contributions due to induced pseudoscalar and other form factors associated with the second class currents come 
into play. In fact all the contributions from the pseudoscalar form factor Ep{Q‘^) and the second class vector form 
factor [Q^) are proportional to while the contribution from the second class axial vector form factor 












11 


is proportional either to ^ or or both. 

To study the lepton mass dependence on r'e(S'e) and induced scattering cross sections in free nucleon as 

well as in nuclear targets, we define A/ = for neutrino/antineutrino induced reaction in and 


e (*^e) 


nuclear targets, where I = i, ii, in, which respectively stands for the cross sections obtained in (i) free 
neutrino/antineutrino-nucleon case, (ii) the Local Fermi Gas Model(LFG) and (Hi) the Local Fermi Gas Model with 
RPA effect (LFG+RPA). 

The results are presented in FiglSJ which show that the differences in the electron and muon production cross 
sections for induced reactions on target are appreciable at low energies E,, < 0.4 GeV. For example, this 

fractional change is about 27(25)% at = 0.2 GeV and reduces to ^ 8(7)% at E^, = 0.4 GeV in the case of free 
nucleon. While in ^‘^G, using LFG it is approximately 40(33)% at E^, = 0.2 GeV and ^ 8 % for both neutrino and 
antineutrino at 0.4 GeV, respectively. However, using RPA effect with LFG, the difference is around 44(42)% and 
^ 9(10)% at El, = 0.2 GeV and 0.4 GeV, respectively. 

While for the case of neutrino(antineutrino) induced process on target using LFG, this fractional change is 
approximately 35(33)% and ^ 7(8)% at E^, = 0.2 GeV and 0.4 GeV, respectively. However, using RPA effect with 
LFG, the difference is around 40(42)% and ^ 9(10)% at Ei, = 0.2 GeV and 0.4 GeV, respectively. 


C. Form factor dependence 

The hadronic current dehned in Eql^ consists of six form factors; three isovector(F% (Q^),i = 1,3) and three axial 
vector{E^{Q^), i = 1, 3) form factors. Among them, {Q^), i = 1,2, are parameterized in terms of Sach’s form fac¬ 
tors G^"(Q^) and G^”(( 5 ^), for which various parameterizations are available in literature leMIl. F^{Q^), which 
arises due to the second class current is generally ignored in calculations. Similarly, the axial current consists of 
three form factors viz. F^{Q^) = Fa{Q'^), = Fp{Q^) and F’^(Q^), among them Fa{Q^) is dominant and is 

parameterized in a dipole form with axial dipole mass(M^). The pseudoscalar form factor Fp{Q^) is given in terms 
of Fa{Q^) and F^{Q^) arises due to second class currents and is generally ignored. In a standard calculation of 
(anti)neutrino nucleon scattering cross sections, form factors associated with first class current are (Q^), F^(Q^), 
Fa{Q^) and Fp{Q^) and their parameterizations are discussed in section|Hl We have used these form factors to calcu¬ 
late cross sections and presented the results in section IHI Al and IHI HI However there are alternate parameterizations 
of the vector form factors [67l - [^ and pseudoscalar form factor and a range of parameter values for the axial dipole 
mass Ma which have also been used in literature to evaluate these cross sections. In the following subsection, we 
give an estimate of the uncertainty in the cross sections associated with the use of alternate parameterizations and/or 
parameter values of these form factors. 


1. Vector form factors 

Vector form factors F^{Q^),i = 1,2 in Eq. (TUI are defined in terms of Dirac and Pauli form factors of proton 
and neutron which in turn are parameterized in terms of experimentally determined Sach’s form factors G^^{Q^) 
and G^"(Q^). For the Sach’s form factors several parameterizations are discussed in literature by various groups like 
Bradford et al. [b^, Budd et al. [b^, Bosted et al. [t^, Alberico et al. [nl , etc. We have studied these parameterizations 
to observe their effect on scattering cross sections for (anti)neutrino induced CCQE processes on nuclear target like 
in local density approximation with RPA effect, i.e. LFG-I-RPA and the results are shown in Fi g El From the 
figure it may be observed that the cross sections obtained using parameterizations of Galster et al. [^, Budd et 
al. [b^, Bosted et al. 0 and Alberico et al. [nl are in a gree ment with each other within 1% which also agrees with 
the recent parameterization discussed by Bradford et al. in the energy region of 0.4 GeV — 0.8 GeV. We have 
also shown the fractional change A/ = (in the inset of FigE]) for vi induced GGQE process using different 

parameterizations and found that the dependence on the choice of parameterizations for F^ (Q^) and F^(Q^) is 
almost negligible. 


2. Axial vector form factor 

The value of axial dipole mass Ma used in Eq. [13] has come recently in debate due to large deviations found 
in the experiments like MiniBooNE (23l - l^ . SciBooNE [2^, K2K [2^, etc. from the world average value 
Earlier measurements for AIa were obtained using vi and vi induced processes on deuterium targets where the 
nuclear effects are expected to play negligible role and the average value is quoted as Ma = 1.02b ± 0.021 GeV 
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FIG. 6: Results of total scatte ring cross section for Ve{ve) induced processes on Ar are shown using different parameterization 
of the vector form factors [gtI-ItI |. Symbol + (Galster et al. [13), dashed dotted line(Bradford et al. [^1. dashed line(Budd et 
al. [S^), dotted line(Bosted et al. and solid line(Alberico et al. [^) are the results obtained using various parameterizations 
for total scattering cross section. Inside the inset A/ = £qj. neutrino induced process using these form factors have 

also been shown. 


Experiment 

Ma {GeV) 

Experiment 

Ma {GeV) 

MINERjzA [^, 

0.99 

SciBooNE [28] 

1.21±0.22 

NOMAD [2^ 

1.05±0.02±0.06 

K2K-SciBar [2^ 

1.144±0.077 

MiniBooNE [23-25] 

1.23±0.20 

K2K-SciFi [29] 

1.20±0.12 

MINOS 

1.19(g'' > 0) 
1.26(g^ > O.SGeV^) 

World Average 

1.026± 0.021 [31J 

1.014±0.014 [3^ 


TABLE II: Recent measurements of the axial dipole mass(Myi). 


while a combined analysis performed by Bodek et al. using Vn.d, and electroproduction data have found 
Ma = 1.014 ± 0.014 GeV. Recent measurements at NOMAD 1^ an d MINER^z A 1^ are close to the world average 
value. On the other hand, the experiments like MiniBooNE |23l - l^ . SciBooNE |28j|. K2K [2^, etc. report a higher 
value of Ma- These experiments were performed with different nuclear targets like as well as some of 

them have used the same nuclear target, for example, NOMAD (2^, MiniBooNE [2^ - 1^ K2K [2^, MINERzzA [2^ 
have used carbon as nuclear target. In Tabl e-HTl w e tabulate the values of axial dipole mass obtained from analyses of 
some recent cross section measurements (2l| - l24 12^.[^. Isol ? ] in the neutrino experiments in the few GeV energy 
range. It is believed that if quasielastic like events which arise due to 2p-2h excitations, meson exchange currents and 
multinucleon correlations are taken into account then the recent experimental results can also be considered to be 
consistent with a smaller value of Ma [TsI - fT^ . However, it may be observed from Table-HIl that even with the same 
nuclear target different values of Ma have been obtained. 

To study the explicit dependence of cross sections on axial dipole mass for neutrino/antineutrino scattering pro¬ 
cesses, we have changed Ma from the base value (taken as the world average value) and obtained the results for 5ma 
defined below in Eq. |M]and Am^i defined in Eq. |3S]by taking the two different values of Ma as 0.9 GeV and 1.2 GeV. 
These results are obtained for the (anti)neutrino induced processes on free nucleon as well as in the LEG with and 
without RPA effect for nuclear target. 

The dependence on axial dipole mass is shown in Figl?! by defining Sma as 


Sma = 


a,, - a,, {Ma = WA) 


ct,,{Ma = WA) 


WA = 1.026 GeV 


( 34 ) 


where I = e ot I = fi. We observe from FigHthat for free nucleon when a modified value of Ma i.e. 
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E^(GeV).E^(GeV)' 


FIG. 7: The dependence of cross section on Ma obtained using Eg. 1341 The results are shown for i^e(t'e) and induced 

processes on free nucleon as well as on Ar target using LFG with and without RPA effect. Solid(dashed) line denotes results 
for the free nucleon case with Ma = 0.9 GeV (1.2 GeV), results obtained using LFG are shown by dashed dotted(double dashed 
dotted) with Ma = 0.9 GeV{1.2 GeV) and results for LFG with RPA effect are shown by dashed double dotted(dotted) with 
Ma = 0.9 GeV{l.2 GeV). 


0.9(1.2) GeV is used instead of world average value of 1.026 GeV then a decrease(increase) of 5 — 15% is obtained 
for Ve/Vfi reactions in the energy range of 0.2 GeV to 0.8 GeV. In the case of Pe/P/x-nucleon reactions this de- 
crease(increase) is about 5 — 10% in the same energy range. When nuclear medium effects are taken into account, for 
example, in the case of Ar nucleus this decrease(increase) remains almost same. Therefore, the uncertainty in the 
(anti)neutrino-nucleus cross sections is the same as in the case of free (anti)neutrino-nucleon scattering processes. 

We now study the sensitivity of the difference in electron and muon production cross sections due to the uncertainty 

= WA) - a^^MA = WA) 
a,JMA = WA) 

Ai(E,)-A2(E,). (35) 

and show the numerical values for Ama for free nucleon and nuclei in Fig. [5] We observe that for free nucleon some 
sensitivity to the difference in the electron and muon production cross sections exists at low energies E^ < 0.4 GeV 
which does not exceed 1%. In nuclear target, when nuclear medium effects are taken into account this sensitivity 
remains almost the same as in the case of free neutrino-nucleon scattering processes. Similar effects are found for the 
case of antineutrino induced reactions. 

From Figs. [7] and El it may be observed that the cross section is very sensitive to the choice of axial dipole mass 
Ma. Therefore, while calculating the charged lepton production cross sections, the value of Ma should be carefully 
taken. 

In the precision era of (anti)neutrino oscillation experiments looking for -o- or o Pe oscillation channels, 
it is important to understand the contribution from the vector part described in Eq. El with E^^ (Q^), iQ^) form 
factors and the axial-vector part with Fa{Q^), Ep{Q^) form factors of the first class current as well as the contribution 


in the choice of Ma- For this we define 

Ai{E,) = 

A2{E,) = 
Ama = 
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FIG. 8: Effect of axial dipole mass on the cross section(from left to right): on free nucleon; LFG, with and without RPA effect 
on target. Here different values of Ma are taken such as 0.9 GeV and 1.2 GeV. The fractional difference(Eq. [SSjl has been 
obtained using the base value of Ma taken as the world average value. Solid(dashed) line denotes results for the neutrino induced 
processes while for antineutrino the results are shown by dashed dotted(double dashed dotted) with Ma = 0.9 GeV(1.2 GeV). 



FIG. 9: Results of the fractional change Spp defined in Eq[3H] as a function of (anti)neutrino energy. The results are shown 
for the Vp induced interaction cross section for the free nucleon case(dashed dotted line), as well as for the nucleons bound in 
'^°Ar(solid line) nuclear target obtained by using LFG with RPA effect. The results corresponding to induced CCQE process 
are shown by dashed double dotted line(free nucleon case) and dashed line("^°Ar target). 


from the second class current with {Q^) and F^{Q^) form factors. Furthermore, the contribution of the terms 
proportional to lepton mass like the contribution from pseudoscalar form factor Fp{Q^) would be different for electron 
and muon channels. To observe their effect on lepton event rates, we have performed calculations using Eqs. [39l l43l 
and szi in which we have studied the individual contribution of the terms in the hadronic current with and without 
Fp(Q2), fY{Q^), F^{Q^) terms. 
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FIG. 10: Results of the fractional change App defined in Eq|^as a function of (anti)neutrino energy. The results are shown for 
the neutrino induced interaction cross section for the free nucleon case(dashed dotted line), as well as for the nucleons bound 
in '*°ylr(solid line) nuclear target obtained by using LFG with RPA effect. The results corresponding to antineutrino induced 
CCQE process are shown by dashed double dotted line(free nucleon case) and dashed line(^°Ar target). 


3. Pseudoscalar form factor 


To study the effect of pseudoscalar form factor Fp{Q‘^) on muon production cross sections, we define 


SppiEp) 


(Pp = 0 ) 

(Fp = 0) 


(36) 


and similar expression for antineutrino is used. For the numerical calculations expression of Fp{Q‘^) given in Eg. 1151 
has been used. The results are presented in Figl9l We find that Spp is more sensitive in the case of induced CCQE 
process than induced process for the free nucleon case as well as for nuclear target. This sensitivity decreases 
with the increase in Vplvp energy and almost vanishes beyond 0.6 GeV. Using the three different expressions of 
the pseudoscalar form factor Fp{Q^) given in Eqs. [13 [TBl and [T71 we have studied the behavior of Fp{Q^) vs 
and found hardly any difference in the dependence and therefore the total scattering cross section for and 
induced processes are also unaffected by the choice of Fp{Q'^). 

We also study the sensitivity of pseudoscalar form factor Fp{Q‘^) to find out the difference in the electron vs muon 
production cross sections that are obtained using EgfTSl For this purpose we define 


Ai(£;,) 


a.,{Fp^Q)-apSFp^Q) 

(PvAFp ^ 0 ) 


(37) 


A 2 (f;,) 


(Fp = 0) - (Fp = 0) 

cTp, {Fp = 0 ) 


(38) 


AFp=A,{Ep)-A2iFp). (39) 

and the results for App are shown in Fisf. 1101 Similar expressions are also used for antineutrino induced processes. 

We have calculated the fractional difference A pp as given in Eq. [53 for free nucleon case as well as for nucleons 
bound in '^'^Ar nuclear target using the LFG with RPA effect. We observe that the inclusion of pseudoscalar form 
factor decreases the fractional change(AFp) by about 3(8)% at ^ 0.2GeU and becomes smaller with the increase 
in energy. When the nuclear medium effects (LFG+RPA) are taken into account in the evaluation of cross sections in 
'^^Ar then this difference increases to 4(15)% at the same energy for neutrino(antineutrino) induced processes. 
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FIG. 11: Results of the fractional change 5pv defined in Eg 1401 as a function of (anti)neutrino energy. The results are shown 
for Vp induced interaction cross section for the free nucleon case(dashed dotted line), as well as for the nucleons bound in 
‘^°Ar(solid line) nuclear target obtained by using LFG. The results corresponding to Vp induced CCQE process are shown by 
dashed double dotted line(free nucleon case) and dashed line('*°ylr target). 


D. Second class currents 




We have also studied the effect of second class current, due to which two additional form factors viz. and 

(Q^) are introduced. 


1. Second class vector current 

The contribution of the second class vector form factor (Q^) to the cross section is always proportional to the 
mass of the lepton so it is quite small in the case of as compared to reactions on free nucleons and nuclei. 

We first study the overall contribution made by the second class vector form factor F^{Q'^) to the cross section in 
the case of muon neutrinos and define 


Spv [Ey) 


^Q)-ay^{FY = Q) 
<yuy{F^ = ^) 


(40) 


Similar expression is used for antineutrino. For the numerical calculations we use Eq. [T8]and the results are shown in 
FigHT] We find that the contribution of F^ {Q^) to the cross section is very small for scattering on free nucleons 
and nuclei. In the case of scattering on nucleons at low energy, the contribution of F^ [Q"^) at Eyjy = 0.2 GeV is 
5% which increases to 7% in when nuclear medium effects are taken into account. 

We now study the sensitivity due to E^ [Q^] in the difference between the electron and muon production cross 
sections for free nucleon and nuclei. Sensitivity that arises in electron and muon production cross sections due to the 
presence of F^[Q"^) is studied by defining 


Ai(f;,) 


A 2 (E,) 

AF3V 


^ ^.y{F^ = ^)-^.AF^ =^) 

truAFA =^) 

= ^i{Ey)-^2{Ey). 


(41) 

(42) 

(43) 


First, we present the results for /S.i{Ey) as a function of neutrino/antineutrino energies and the results are shown in 
Fig. [T^ These results are presented for the free nucleon case as well as bound nucleons in Ar and the cross sections 
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FIG. 12: The difference of fractional changes Ai defined in Eas llll for the free nucleon case(solid line), and in for 

neutrino(Left panel) and antineutrino(Right panel) induced processes using LFG(dashed line). 


are obtained using the expression of the hadronic current with second class vector current. We must point out that 
the contribution from the second class axial current is switched off. We find that Ai{E^) is sensitive to the flavor of 
neutrinos especially at low energies [Evjv < 0.3GeV) which is mainly due to threshold effect. When we perform 
calculations on nuclear targets like using LFG, the fractional difference Ai{E^) changes from free nucleon case. 
However, it is not sensitive to the choice of nuclear model. 

In Fig. [131 ■we present the results for Apv using Eq. |43]and obtained it for the free nucleon case as well as for 
and '^^Ar nuclear targets. We find that the effects are energy dependent and are more pronounced at low energies. 
From the figure it may be noticed that the fractional change is the same for both and nuclei. For example, for 
the case of neutrino, at Ei, = 0.2 GeV, Apv is ^ 1% for free nucleon as well as in '^^Ar evaluated using LFG and the 
difference becomes almost negligible beyond E^, = 0.5 GeV. Similarly, for the case of antineutrino, at E^ = 0.2 GeV, 
Apv is ~ 4% for the free nucleon case as well as in '^^Ar evaluated using LFG, which becomes ~ 1% at Ei, = 0.5 GeV. 

In general there is a large uncertainty associated with the determination of (Q^) and F^{Q^) form factors. 
We have also studied uncertainty due to various parameterizations of form factor {Q^)- Some of the alternative 
parameterizations of the form factor F^ {Q^) are given in Eqs|THl[THl which have been used for numerical calculations. 
The results are shown in FigjTd] We find that for neutrino induced process on free nucleon target, the difference in 
the results for AF^ obtained by using two different forms of F^ {Q"^) (Form I using Eq. |T51 Form II using Eg. IT^ is 
very small. For example, this difference is ^ 1% at low energies(^ 0.2 GeV) which almost vanishes with the increase 
in energy. In the case of antineutrino induced reaction on free nucleon, this difference is around 3% at low energies 
which gradually vanishes with the increase in energy. 


2. Second class axial current 


The axial vector form factor associated with second class current(F^(Q^)) also contributes to the cross section in 
addition to the second class vector form factor(F;^(Q^)). To observe the effect of F^{Q‘^) on induced cross 

section we define fractional difference 


5FA{Ey) 


^..(^ 3 ^= 0 ) 


(44) 


and similar expression for antineutrino is used. We show the numerical results in Fig llSI These results are presented 
for the free nucleon case as well as for bound nucleons in '^^Ar and the cross sections are obtained using the expression 
of the hadronic current with second class axial vector current. This is to point out that the contribution from the 
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FIG. 13: The difference of fractional changes defined in Eg 1431 for the free nucleon case(neutrino results shown by dashed- 
dotted line and antineutrino results by dashed-double dotted line) as well as for (circle for neutrino and triangle up for 
antineutrino) and Ar (solid line for neutrino and dashed line for antineutrino) nuclear targets obtained by using LFG. 



FIG. 14: Variation of as given in Eqs. |43]as a function of (anti)neutrino energies are shown for different forms of (Q^) 
used in second class currents. We have performed the calculation for (anti)neutrino induced scattering processes on free nucleon 
target. The results are obtained for Form I and Form II using Eos. 1181 and 1191 respectively for Ma = 1.026 GeV. The results 
obtained by using Form I for neutrino(antineutrino) are shown by dashed-dotted (dashed-double dotted)line and with Form II 
are shown by solid(dashed)line, respectively. 


second class vector current is switched off. We find that 5pA is hardly sensitive to the presence of in and 

scattering cross sections from free nucleon and nuclear targets. We also study the sensitivity of the electron and 
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FIG. 15: Results of the fractional change SpA defined in Eg 1441 as a function of (anti)neutrino energy. The results are shown 
for the induced interaction cross section for the free nucleon case(dashed dotted line), as well as for the nucleons bound in 
‘^°Ar(solid line) obtained by using LFG. The results corresponding to induced CCQE process are shown by dashed double 
dotted line(free nucleon case) and dashed line('*°ylr target). 


muon production cross sections to F^{Q^) for free nucleon and nuclei, by defining 




^2{Eu) 


^uAEi = 0) 

- d.2{E,) 


(45) 

(46) 

(47) 


and the numerical results are shown in Fig ll61 for We find the results to be similar in nature as found in the 

case of currents with second class vector form factor. In Fig.[T71 we present the results for /S.pA. In this case also we 
find the sensitivity to be smaller than observed in the case of A^v. 

When we compare our present results for the difference in the electron and muon production cross sections on 
free nucleon target with the results of Day and McFarland [l^, we find that our results for the contribution of the 
F^{Q^) in the case of antineutrino reactions and the results for the contribution of F^{Q'^) for neutrino reactions 
agree qualitatively with their results. This is not so in the case of F^[Q"^) for neutrino reactions and ^ 3 ^( 15 ^) for 
antineutrino reactions. This may be due the different expressions used for the contribution of the interference terms 
between first and second class currents. Our expressions agree with the general expressions given by Pais [Slj, Kuzmin 
et al. [13 but not with Eq.(3.18) of Llewellyn Smith [4^ which has been used in Ref. [ 2 ^. 


E. Radiative corrections 


Radiative corrections are potential source of difference between electron and muon production cross sections in 
(anti)neutrino reactions due to their logarithmic dependence on the lepton mass through terms like log(-^), where 
E* is some energy scale in the reaction. The radiative corrections in the charged current quasielastic neutrino-nucleon 
reactions relevant for the present oscillation experiments in the energy region of few GeV have been recently calculated 
by Bodek [s^. Day and McFarland [ 2 ^ and Graczvc k f^ . Bodek [83| and Day and McFarland make use of 
leading log approximation given by De Rujula et al. [83 to calculate the contribution of soft photon emission by 
the lepton leg bremsstrahlung diagram which gives major contribution to the radiative corrections depending on the 
lepton mass mi. On the other hand, Graczyck 8^ includes the contribution of other diagrams like two boson exchange 
involving W and 7 , propagator correction in addition to the soft photon bremsstrahlung. We have used the results 
of De Rujula et al. (85l | for the radiative corrections in the neutrino nucleon scattering and study the effect of nuclear 
medium on these radiative corrections. In the work of De Rujula et al. [85l| . the modified cross section including 














20 



FIG. 16: The difference of fractional changes Ai defined in Eos 145 1 for the free nucleon case(solid line), and in for 

neutrino(Left panel) and antineutrino(Right panel) induced processes using LFG(dashed line). 



FIG. 17: The difference of fractional changes defined in Eg 1471 for the free nucleon case(neutrino results shown by dashed- 
dotted line and antineutrino results by dashed-double dotted line) as well as for (circle for neutrino and triangle up for 
antineutrino) and (solid line for neutrino and dashed line for antineutrino) obtained by using LEG. 


radiative corrections is given by 


da 

dEidfli 


da free O: f 1 4- .2: / 1 da free f./ \ da free \ 


( 48 ) 


















21 



E,(GeV) 

FIG. 18: The effect of radiative corrections on fractional difference Arc defined in Eq. [5T]for (anti)neutrino induced processes on 
free nucleon as well as on target using LFG with and without RPA effect. For the neutrino(antineutrino) induced processes 
on free nucleon the results are shown by dashed dotted line(plus), results for calculations using LFG for neutrino(antineutrino) 
are shown by dashed line (triangle up) and results for calculations using LFG with RPA effect for neutrino(antineutrino) are 
shown by solid line (circle). 


where a free is the (anti)neutrino induced cross section obtained without radiative effects, E* is the lepton energy in 
the center of mass of neutrino nucleon system and Zmin is given by 


^min 


with E* 


AEf 
2ME^ 
s + mf — 

Vi 


(49) 

(50) 


where s = {k + p)'^ is square of total energy in the center of mass system, k and p are the four momenta of incoming 
neutrino and target nucleon. 

To show the effect of radiative corrections on the lepton event rates, we have obtained total scattering cross sections 
for Ve and induced reactions on free and bound nucleon, with and without radiative corrections and define 


^i{E,) 

^2{E,) 

^RC 


a,^{RC) - a,ARC) 

(Ju, (RC) 

a,^{NR)-a,SNR) 

<y.ANR) 

^l(V) ~ ^2(V) 


(51) 


where a^^{RC), (/ = e,/i) represents the cross sections obtained by taking radiative corrections into account and 
a^i{NR), {I = e,p) represents the cross sections without radiative corrections. A similar definition has been used 
for antineutrino induced process. The results for Arc are presented in Fig. [181 These results are shown for the 
free nucleon target as well as for per interacting nucleon in '^^Ar, for which the total cross sections are calculated 
in LFG with and without RPA effect. We find that the difference in electron and muon production cross section 
Arc due to radiative corrections is quite small in case of neutrino induced reactions as compared to antineutrino 
reactions. Furthermore, the effect of nuclear medium is to further decrease Arc and the reduction is larger in case 
of antineutrino reactions as compared to neutrino reactions. For example, for the neutrino induced process on free 
nucleon the effect is around 2% at E^, = 0.5 GeV while for antineutrino the effect is around 5% at E^, = 0.5 GeV. 
When we performed the calculations using LFG with and without RPA effect, the effect becomes less than 1% in the 
case of neutrinos. For antineutrino induced process, using LFG without RPA effect, it is around 3%, which becomes 
~ 2% when RPA effect is included. 
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IV. SUMMARY AND CONCLUSIONS 

In the present study we observe the following: 

1 - induced cross sections in free nucleon as well as nucleons bound in nuclear targets are more suppressed due 

to threshold effects at low energies than Ve{ve) induced reaction cross sections. Moreover, when cross sections 
are evaluated in nuclear targets there is a further reduction in the cross sections due to nuclear medium effects. 
This reduction is energy dependent. For example, when calculations are performed for nucleus in LFG, the 
cross section for Ve{ve) induced scattering is reduced by ^ 30(45)% at Ey = 0.3 GeV from free nucleon case 
while for induced scattering the reduction is ~ 32(46)% at Ey = 0.3 GeV. Inclusion of RPA correlation 

with LFG, further reduces the cross section from J^e(j^e) induced scattering by ^ 37(25)% at Ey = 0.3 GeV, 
which for v^{p^) induced scattering is ^ 40(27)% at Ey = 0.3 GeV. This results in a larger difference in electron 
and muon production cross sections for the case of nuclear targets as compared to the free nucleon target. 

2. At low energies of Ey/y < 0.5 GeV there is appreciable nuclear model dependence on (anti)neutrino-nucleus 
cross sections for both flavors of neutrino (antineutrino). The suppression due to nuclear medium effects is 
larger in the Local Fermi Gas Model(LFG) as compared to the Fermi gas model of Llewellyn Smith [d^. The 
suppression in the Fermi gas models of Smith and Moniz [d^ and Gaisser and O’Gonnell [4^ are larger than 
LFG. When RPA effect is included in LFG, the suppression is largest. 

3. The suppression due to nuclear medium effects is larger in the case of antineutrino as compared to the cross 
sections obtained for neutrino induced processes. 

4. For a given set of parameters which determine the form factors and other coupling constants the percentage 
difference in electron and muon production cross sections is more for nuclear targets than for the free nucleon 
target. This difference decreases with neutrino/antineutrino energy. Also this difference increases with the 
increase in mass number. 

5. The percentage difference in electron and muon production cross sections due to uncertainties in axial dipole 
mass is more in the case of nuclear targets as compared to free nucleon target. The difference increases with 
the increase in mass number. 

6 . The fractional difference in the cross sections due to the presence of pseudoscalar form factor is more in the case 
of induced GGQE process than induced process for the free nucleon case as well as in nuclear targets. 
Qualitatively at low neutrino energies there is a small difference in the results obtained for the free nucleon target 
and nucleons bound in nucleus. This difference vanishes with the increase in energy. In the case of antineutrino 
induced reaction the difference is slightly larger than found in the case of neutrino and this difference does not 
vanish with the increase in energy. The difference is almost independent of the choice of nuclear target. 

The contribution of pseudoscalar form factor in nucleon scattering is about 3(9)% at low energies of 

Ey = 0.2 GeV and becomes ^ 1(3)% at Ey = 0.4 GeV in the case of free nucleon. When nuclear medium effects 
are taken into account this contribution increases. 

7. The inclusion of second class vector current results in an increase in the total scattering cross section if present 

experimental limits of the second class form factor (Q^) is used for r'e(j^e) as well as v^{p^) induced processes. 
With the inclusion of (Q^), induced scattering cross section increases about 1(4)% for free nucleon 

case and 1(8)% when we include RPA with LFG at Ey = 0.4 GeV. For the Ve{^e) induced processes this effect 
is smaller than in comparison to the induced processes. The difference in electron and muon production 

cross sections increases at low energies for nuclear targets as compared to the free nucleon target. This difference 
almost vanishes with the increase in neutrino energy. In the case of antineutrino induced GGQE process this 
difference is slightly more as compared to the neutrino case. However, this difference is almost independent of 
the choice of nuclear target. 

The effect of including second class axial vector form factor F^{Q^) (consistent with present experimental limits) 
is qualitatively similar to the effect of including second class vector form factor but quantitatively quite small 
as compared to the effect of F^ (Q^) summarized above. 

8 . The effect of radiative corrections being proportional to log{—) affects the r'e(Pe) scattering cross section more 
than v^{v^) scattering cross sections when a corresponding charged lepton is in the final state. This gives a 
difference in z^e(Pe) and Vfi{vfi) scattering cross sections which is almost independent with energy in the case of 
neutrino induced process while it increases slightly with energy in the case of antineutrino induced process both 
for the free nucleon as well as the bound nucleons. We find this difference to decrease in the presence of nuclear 
medium effects and this decrease is more in the case of antineutrino as compared to the neutrino reactions. 
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Appendix A: Hadronic tensor J'"' 

1. Relativistic expression for the hadronic tensor J^’' 

We have followed the prescription of Nieves et al. [i^ while taking RPA correlations into account. The difference is 
that we have also included the contribution coming due to the second class currents. The following is the expression 
for the hadronic tensor J^'^: 


Ji^'' = 4 


+ q^p" + ^P^p" + q^g^’'/^) - {f^? \ + (i. 


\ 


AM^ 




jA + 2p»p-)| + {F^ f EmV'- I - 1 I + + pV + 2p“p’' 

2 


+ F^FY {q^g^^'' - qf^q^) - 2 iiF^F a + F^FA)e^‘'‘^^pp,qp - 2 FAFpq^q'^ 
- F^Fp^ip^^q'' + q'^p'' + q^q’^) - 2 F,^Fa [p^'' + gV + q^ql 
+ 2FYF^ ip^q’^ + q^p'' + q>^q'^) + F^F^^{p^q’^ + q^p‘' + q'^q’') 


(Al) 


2. Component form of incorporating RPA in the lowest order 

Here the three momentum transfer q is taken along z-axis and RPA is applied in the leading terms. The different 
components of the hadronic tensors with RPA effect are 


7 OO _ A 
•^RPA — ^ 


.iFir (c^E^ip) P £ + ,»E(p)) - £(riT f <A . + + 


q" 


M2 


+ 2{FA f ^ _^(Fp)2|^ + 2{Flf{qy (^1 - 

- {F^f^ ~ + 2^1"^^3"" {{q°f + 2A(p)g°) 

^ - ‘^FpFA{q°r - ^FaF^ {E{p)q° + {q°f) 


+ F 


- 2FpF^ 


,A q 


{E{p)q^ + {qy) 


(A2) 
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tOz 


RPA 


= 4 


(F^ f {C^Eip)i2p, + |g I) + gV) - 


f E{p) 2p^ + |g I ^ g°|g| q°{2p^ + |g I) 


M 


2 M 2 


+ (Fa)^ {Ci^E{p){2p, + |g1) + gV) - C^{Fp)\°\c^^ + 2{F^)\°\q] (^1 - 
- (^3^)"^(I91 + 2pz)(g° + 2E{p)) - F^F^q^Wl + 2F^F^{E{m + gV + q°m 
+ FYFY7^{E{iY\^ + q°Pz + q°\q\) - 2FpFAq°\^ - 2 FaFY{E{p)\^ + q°Pz + g°M) 


2 3 2 M 2 ’ 

r,2 


- i^pF. 
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M 2 
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4 y 2 
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2M 
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0^2 
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2 

lAI I I l?V ipV 'i 


- + b1)' - + ^FYfY{2pM + f) + fYfY^{2pM + f) 


- FpFY^{2pz\q\ + if) - 2FAFpf - 2 FaFY{2pM + f) 


rxx _ A 
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2{FY)HpI- 
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pI - gV4 
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- 2{FY)^^-^fFYFY 


M 2 
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(A3) 


(A4) 

(AS) 

(A 6 ) 


Cn = 


1 - ' 2 o/'(p)C^Ar(g, fcp)P’ 


Ct = 


l-F(g,fcp)Ft(g)| 2 ^ 


Cl = 


l-C/(g,fcp)Cz(g )|2 


where Vi and Vt are the longitudinal and transverse part of the nucleon-nucleon potential calculated with tt and p 
exchanges and are given by 


Vl{q) = ^ 


Mq) = 



t 

47r 


= 0.08, A^ = 1.2 GeV, = 0.14 GeV 


+ g' } , Gp = 2,Ap = 2.5 GeV, rup = 0.77 GeV 


(A7) 


g' is the Landau-Migdal parameter taken to be 0.7 which has been used quite successfully to explain many electromag¬ 
netic and weak processes in nuclei. ?7(g, fcp) = U]s[{q,kp) + U/^{q,kp) is the Lindhard function for the particle-hole 
excitation and U/^{q,kp) is the Lindhard function for the delta-hole excitation. The details are given in Ref. [i^. [ 
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